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ABSTRACT 

To move one step forward toward a Galactic distribution of planets, we present the first planet sensitivity 
analysis for microlensing events with simultaneous observations from space and the ground. We present this 
analysis for two such events, OGLE-2014-BLG-0939 and OGLE-2014-BLG-0124, which both show substan¬ 
tial planet sensitivity even though neither of them reached high magnification. This suggests that an ensemble 
of low to moderate magnification events can also yield significant planet sensitivity and therefore probability to 
detect planets. The implications of our results to the ongoing and future space-based microlensing experiments 
to measure the Galactic distribution of planets are discussed. 

Subject headings: gravitational lensing: micro — stars: planet 


1. INTRODUCTION 

Not relying on the light from the target system, microlens¬ 
ing is in principle sensitive to planets at various line-of-sight 
distances, suggesting that a sample of microlensing planets 
will be able to tell us the Galactic distribution of planets. 
However, a problem of the standard microlensing technique 
is that the mass Ml and distance of the lens system, to¬ 
gether with the lens-source relative proper motion /x re 1 , are 
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buried within a single observable, the event timescale t e, 
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Although a few methods have been used to resolve this de¬ 
generacy, they are either ineffective, in the sense that they 
can only be applied in very rare cases (Gould et al. 2009; 
Yee et al. 2009; Gould & Yee 2013), or strongly biased to¬ 
ward nearby lenses or long-timescale events (Gould 1992). 
Therefore, previous statistical studies based on microlensing 
planets had to assume some typical values of the lens system. 
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for example, Z>l = 4 kpc and Ml = 0.3M 0 (e.g., Gould et al. 
2010; Clanton & Gaudi 2014). This is appropriate for com¬ 
plementing the demographics of planets by combining with 
other detection techniques (Gaudi 2012), but it has prevented 
microlensing from demonstrating its unique power: deriving 
the distribution of planets at various Galactic distances. 

This situation has been changing with the emergence of 
space-based microlensing programs. As has long been re¬ 
alized, combining observations from the ground and at least 
one satellite that is well separated (~ AU) from Earth is the 
only effective way to measure the mass and distance of poten¬ 
tially all microlenses (Refsdal 1966; Gould 1994), because the 
measurable parameter from such an experiment, the microlens 
parallax vector 7 Te, with amplitude and direction defined by 


Ttrel _ 7TK _ /Tel 
0e ’ 7Te Hiel 


( 1 ) 


relates to Ml and Z>l by 


$E 1 _ 1 + 7T E (9e 
ktte ’ D l D s AU 


( 2 ) 


With observations from Earth and only one satellite, there are 
four allowed choices of 7Tp. This four-fold degeneracy arises 
from the fact that the projected position of the source onto the 
sky as seen from Earth and the satellite can pass by the lens 
from the same side or different sides. Given that the majority 
of Galactic microlensing events have well measured source 
distance Ds (~ 8 kpc), Equation(2) suggests that dermining 
M l and D] requires measurements of both 7 t e and Op. For 
events with detected planets, 6e can usually be measured via 
finite source effects (Yoo et al. 2004). However, measuring 
the Galactic distribution of planets requires a determination 
of the distance distribution of all the lenses being probed, 
not just those with planets. Or, more precisely, the distance 
distribution of the planet sensitivity of the events containing 
these lenses. For single-lens events, only a small fraction 
can yield definite distance measurements even though space- 
based parallax can be measured (Zhu et al. 2015b). Fortu¬ 
nately, Calchi Novati et al. (2015) showed that individual dis¬ 
tances can, in the great majority of cases, be statistically in¬ 
ferred from a combination of 7 Te measurements and kine¬ 
matic priors derived from a Galactic model. The 2014 pilot 
program using Spitzer as a microlens parallax satellite has 
demonstrated the feasibility of such a strategy to derive the 
Galactic distribution of planets (Calchi Novati et al. 2015), a 
key ingredient of which is measuring parallaxes to constrain 
or measure the lens (host star) distances. 

To conduct statistical studies with these microlensing 
events, we will also need a good understanding of the de¬ 
tection efficiency of planets in each individual event. The 
estimation of detection efficiency, or planet sensitivity, via 
microlensing has been studied for a long time, both ana¬ 
lytically and computationally. Soon after microlensing by 
planetary systems was proposed by Mao & Paczynski (1991), 
Gould & Loeb (1992) estimated that about 20% microlens¬ 
ing events would show planetary anomalies if all lenses 
are solar-like systems. A planet may be detected in a mi¬ 
crolensing event if the source (or at least some portion of it, 
Ingrosso et al. 2009) passes over or near a caustic induced by 
the planet. Hence, the position of the source at a particular 
time determines what planets the event is sensitive to at that 
moment. The total planet sensitivity of an event is an integra¬ 
tion of the sensitivity over the light curve. The light curve may 


also be broken into segments each with its own planet sensi¬ 
tivity. The methodology to compute the planet sensitivity was 
first proposed and used in Rhie et al. (2000); Gaudi & Sackett 
(2000) and Albrow et al. (2000). 

After the search networks were established, microlensing 
began giving meaningful constraints on the planet occurrence 
rate. For example, Gaudi et al. (2002) put a 33% upper 
limit on the occurrence rate of Jupiter-mass planets around 
Bulge M dwarfs after analyzing 43 intensively monitored 
events; Gould et al. (2010) for the first time presented the 
planet frequency beyond the “snow” line, using an ensemble 
of 13 extremely high-magnification (A > 200) events. The 
planet frequency from microlensing has also been studied by 
Sumi et al. (2010) and Cass an et al. (2012). 

The high-magnification events in the Gould et al. (2010) 
sample are the most efficient for measuring planet frequency 
because they are each very sensitive to planets, achieving 
up to 100% detection efficiency for the largest mass ratios 
(Griest & Safizadeh 1998). However, such events are ex¬ 
tremely rare and difficult to predict in real-time (which is nec¬ 
essary for both ground-based characterization and scheduling 
Spitzer-type observations). Given these difficulties, it is nat¬ 
ural to turn to the much larger sample of low-magnification 
events to measure the planet frequency. While each one is far 
less sensitive to planets, their number makes up for this defi¬ 
ciency, allowing for a measurement of the planet frequency as 
in Gaudi et al. (2002) or Cassan et al. (2012). 

The difficulty with low-magnification events is that the ar¬ 
guments used by Gould et al. (2010) to establish an unbiased 
sample do not apply. However, Calchi Novati et al. (2015) 
suggested that any event can be included in the sample to mea¬ 
sure the Galactic distribution of planets as long as they are 
selected for space-based parallax observations without refer¬ 
ence to the presence or absence of planets. This idea has been 
codified for Spitzer and for future space-based microlensing 
experiments in Yee et al. (2015b). In brief, events may be 
chosen either objectively or subjectively. An event is ob¬ 
jectively selected if it meets some pre-determined criteria, in 
which case, the planet sensitivity from the entire event may be 
considered. For a subjectively selected event, only the planet 
sensitivity from the portion of the light curve observed after 
the selection is made may be considered. This is necessary 
to prevent prior knowledge about the event from influencing 
the planet frequency measurement (see Yee et al. 2015b, for a 
more detailed discussion). 

In this work, we present the planet sensitivity analysis of 
two microlensing events monitored in the 2014 pilot program. 
The purpose is to illustrate the various issues that arise in 
this kind of analysis. We give an overview of our compu¬ 
tational method in Section 2, and present the results, i.e., the 
planet sensitivities of two events OGFE-2014-BFG-0939 and 
OGFE-2014-BFG-0124, in Section 3. Some discussion on 
the implications of these results is given in Section 4. 

2. METHODS OVERVIEW 

The procedure for computing the planet sensitivity for mi¬ 
crolensing events with ground-based and space-based obser¬ 
vations has been outlined in Section 4 of Yee et al. (2015b). 
We further improve the methodology below by filling in more 
details about the computation. 

We first measure the planet sensitivity as a function of two 
parameters, S(q,s), where q is the planet-to-star mass ratio 
and .s' the projected separation normalized to the Einstein ra¬ 
dius 0e . The other parameters required to describe a planetary 
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Fig. 1.— Example light curves (upper panels) and corresponding trajectories seen from Earth and Spitzer (lower panel) generated for the event OGLE- 
2014-BLG-0124. The planet-to-star mass ratio q = 7 x 10 -4 and the projected separation s = 0.94 of the detected planet OGLE-2014-BLG-0124Lb are used 
(Udalski et al. 2015). In the upper panel, the colors encode which team took the observation at that time, and the typical uncertainties of observations taken by 
different teams are indicated by the error bars to the lower right, In the lower panel, the planet and its host are placed at (0.94,0) and (0,0), respectively, and the 
caustic arising from this lens configuration is shown in the black curve. The arrows indicate the directions of source relative motion for each a value. In fact, the 
trajectories seen from Earth and Spitzer are fixed by the measured 7 Te and uq and have directions almost due east. While in reality, the orientation of the source 
trajectory is fixed on the sky and the light curves shown represent different orientations of the caustics, we choose to plot multiple sets of trajectories rather than 
multiple caustics in order to avoid severe overlapping among caustics. 


event include: the time of maximum magnification to, event 
impact parameter no, the Einstein timescale the scaled 
source size p, and the planet-star axis orientation a (relative 
to the source-lens trajectory). 

We adopt the approach that was first proposed by Rhie et al. 
(2000) to compute S(q,s). For each set of ( q,s ), we generate 
300 planetary light curves that vary in a but have other pa¬ 
rameters the same. Besides ( q,s ), we adopt (fo, woTe) from the 
best fit of the single-lens/planetary event. We adopt the value 
of p if it is measured via the finite-source effect, which is typ¬ 
ical in planetary events (such as OGLE-2014-BLG-0124) or 
extremely high-magnification single-lens events. Otherwise, 
we choose p as Yee et al. (2015b) suggested. That is, combin¬ 
ing f E , 0*, and a reasonable choice of // re | to determine p by 
p = (9*/0e = $*/(fii-eif E ), as in the case of OGLE-2014-BLG- 
0939. 5 

5 Note that while a fairly precise estimate of p is often required to properly 


Each planetary light curve is generated by creating fake data 
points at the times when the real measurements were taken, 
with values equal to those predicted by the model and error 
bars that are the same in magnitudes as those of the real data 
points. To maximize the efficiency of the planetary light curve 
computation, we use the point-source, quadrupole, and hex- 
adecapole (Pejcha & Heyrovsky 2009; Gould 2008) approxi¬ 
mations when the source is approaching but still reasonably 
far (> 2 p) from the caustics. For epochs that are near or on 
crossing caustics, we use contour integration, in which the 
limb-darkening effect is accommodated by using 10 annuli 
(Gould & Gaucherel 1997; Dominik 1998). This contour in¬ 
tegration may fail under specific lens-source configurations 
(for example, if the lens sits on a sharp cusp, the contour 
integration may fail to identify some tiny images), in which 
case the more time-consuming inverse ray-shooting is used 

model real planetary microlensing events, only a rough estimate is needed to 
estimate the planet sensitivity of point-lens events (Gaudi et al. 2002). 
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(Dong et al. 2006). 6 We then find the best-fit single-lens 
model of each of these fake light curves using the downhill 
simplex algorithm, and quantify the deviation between the 
best-fit model and the fake data by Xsl- 7 In searching for 
the best-fit single-lens model, any unphysical model that has 
severely negative blending (h H < -0.2, i.e., an I < 19.75 “anti¬ 
star”, Smith et al. 2007) is automatically rejected. Here F = 1 
is defined to correspond to I = 18. With the absence of obser¬ 
vational systematics and statistical fluctuations in these fake 
light curves, any x| L would be due to the presence of the in¬ 
jected planet. We choose x| L = 200 as the threshold for a 
“detection”, which has been shown to be reasonable by vari¬ 
ous simulations (e.g., Zhu et al. 2014; Henderson et al. 2014). 
Then the fraction of all the a at fixed (q. s) for which the planet 
is detectable is said to be the sensitivity S{q,s). 

In principle, one might argue that the criterion we set above 
for claiming a “detection” is only reasonable in simulated 
events but not enough for real planet detections, because it 
is possible that most of the \ 2 ma y come from a single or two 
data points. This is especially true for small planets in events 
with low cadence, such as OGLE-2014-BLG-0939 (Yee et al. 
2015a). However, in the presence of various systematics in 
real microlensing observations, such an anomaly cannot be 
securely claimed as a detection. Thus, it may seem that other 
criteria are required to take this issue into account. For exam¬ 
ple, Shvartzvald &Maoz (2012) requires at least three con¬ 
secutive data points each with a 3 a deviation from the best-fit 
single-lens model to simulate the anomaly detection process 
in survey mode (e.g.,the OGLE Early-Early Warning System 
anomaly detector, Udalski 2003). The problem of such a crite¬ 
rion is that it does not take into account the presence of poten¬ 
tial follow-up observations. Once a small anomaly is found 
in the survey data, the follow-up teams will obtain more in¬ 
tensive observations to confirm the nature of this anomaly, 
but because these follow-up observations are triggered by the 
anomaly, they cannot be included in estimating the planet sen¬ 
sitivity. Therefore, we think our single criterion based on y 2 
remains reasonable in a more realistic situation. 

For each event, we first compute S(q,s ) in a grid with ten 
q values equally spaced in log r/ from 1CT 5 to 10 -2 and ten ,v 
values equally spaced in logs from 0.1 to 10. Based on this 
preliminary S{q,s ) map, a finer grid is constructed and the 
final S(q,s) map is computed. In the end, the planet sensi¬ 
tivity as a function of q is obtained by marginalizing over all 
s values, and the integrated planet sensitivity S lo , is given by 
further marginalizing over all q values, 

S(q)= JS(q,s)dlogs; S tot = J S(q)£(q)d\ogq . 

The first step assumes a flat distribution in logs (i.e., quasi 
Opik’s law), and the second step adopts some planet mass 
function £,(q). 

As has been emphasized in Yee et al. (2015b), only obser¬ 
vations that were carried out without reference to the presence 
or absence of a planet can be included when S(q,s ) is being 

6 Penny (2014) developed a new algorithm to speed up the computation of 
light curves with extremely low-mass planets (q < 10 5 ). Because of the lack 
of such planets in our current analysis, we end up not using this algorithm. 

7 The other approach to compute S(q,s ), first used in Gaudi & Sacked 

(2000), requires fitting planetary models to the real data, and thus is extremely 
time-consuming, although it has the advantage of simultaneously searching 
for all planets that may be lurking in the data down to the adopted threshold 
(Gaudi et al. 2002) 


TABLE 1 

Microlensing parameters of OGLE-2014-BLG-0939 and 
OGLE-2014-BLG-0124 USED IN THE PLANET SENSITIVITY 
CALCULATION. PARAMETERS OF OGLE-2014-BLG-0939 ARE TAKEN 
FROM Yee ET AL. (2015a), AND PARAMETERS OF 
OGLE-2014-BLG-0124 ARE FROM UDALSKI ET AL. (2015). 


Parameters 

Unit 

OGLb-2014-BLG-0939 
(mo solution) 

OGLE-2014-BLG-0124 
(mo > 0 solution) 

to - 6800 

day 

'36.20 

3076 

uo 

— 

-0.913 

0.1749 

tE 

day 

22.99 

152.1 

7TE,N 

— 

0.220 

-0.0055 

te.e 

— 

0.238 

0.1461 

P 

— 

0.012 

0.00125 

9 

— 

— 

0.00694 

s 

— 

— 

0.9443 


calculated. The treatment of this constraint varies for differ¬ 
ent events. We return to this point with individual examples 
in Section 3. 

We apply our method to the first two events that have simul¬ 
taneous observations taken from space and ground, OGLE- 
2014-BLG-0939 and OGLE-2014-BLG-0124. 8 For a sim¬ 
ple demonstration of our method, we show in Figure 1 some 
fake light curves generated for the planetary event OGLE- 
2014-BLG-0124. The injected planet has q = 7 x 1CT 4 and 
s = 0.94, both taken as the detected values. Readers can find 
the real OGLE-2014-BLG-0124 event light curve in Figure 2 
of Udalski et al. (2015). 

3. RESULTS 

Below we present the planet sensitivity analysis of two 
events, OGLE-2014-BLG-0939 and OGLE-2014-BLG-0124, 
both from the 2014 Spitzer microlens parallax program. Al¬ 
though follow-up observations in most cases are crucial to 
characterize the properties of microlensing planets, they are 
usually taken as responses to the anomaly in the data taken 
by survey teams in the case of low or moderate magnification 
events. Since only those observations that are taken without 
reference to the presence or absence of planets can be used to 
derive the planet sensitivity, we only include survey data in 
the present analysis. 

The microlensing parameters of these two events, used for 
planet sensitivity computations, are given in Table 1. 

3.1. Single-lens event OGLE-2014-BLG-0939 

OGLE-2014-BLG-0939 was the first single-lens event 
with space-based microlens parallax measurement (Yee et al. 
2015a). The four-fold degeneracy, arising from the fact that 
the projected position of the source onto the sky as seen from 
Earth and the satellite can pass by the lens from the same side 
or different sides (Refsdal 1966; Gould 1994), was effectively 
broken given the measurement of the source proper motion, 
yielding a star with mass M = 0.23 Mq at distance I)\_ = 3.1 
kpc. However, since the four-fold degeneracy is a generic 
feature of most single-lens events, we also consider below the 
case in which it would not have been broken. 

The event OGLE-2014-BLG-0939 was only observed 
by the Optical Gravitational Lensing Experiment (OGLE, 
Udalski 2003). It lies in OGLE field BLG630, implying that 

8 Although OGLE-2005-SMC-001 is the first microlensing event with a 
space-based microlens parallax measurement (Dong et al. 2007), it is not in¬ 
cluded here, because its Spitzer light curve has too sparse observations to 
allow a planet sensitivity analysis, and because the lens system is a stellar 
binary. 
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Fig. 2.— Planet sensitivity curves of OGLE-2014-BLG-0939. The solid, dashed, dot-dashed, and dotted curves represent the sensitivity curves from 
“OGLE+Spitzer,color”, “OGLE+Spitzer’, “OGLE only”, and “OGLE+Spitzer,alert”, respectively. In all these cases the ) solution, i.e., the correct solu¬ 
tion, has been used. Panels (a), (b) and (c) show the “OGLE+Spitzer,color” map and one of the other three sensitivity maps, with the colors representing the 
curves with different sensitivities in S(q,s), and the horizontal dashed lines indicating the masses of three Solar system planets. The projected separation ,s = 1 
corresponds to a physical separation of 1.9 AU. Panel (d) shows the four S(q) sensitivity curves in a single plot. 


it was observed at relatively low cadence, roughly once per 
two nights. We include OGLE data taken from HJD' = HJD 
-2450000 = 6000 until the end of the 2014 season (HJD' ps 
6941). These boundaries are somewhat arbitrary, but allow 
for a substantial baseline while keeping the number of data 
points manageable. After the removal of isolated outliers, we 
find in total 248 observations. This event also received in total 
31 Spitzer observations during the interval HJD' = 6814.1 to 
6845.7, with close to once per day cadence. Readers can find 
more details about these observations in Yee et al. (2015a). 

Although event OGLE-2014-BLG-0939 was not expected 
to show substantial planet sensitivity because of the low ob¬ 
servational cadence, it is chosen for the present analysis for 
other reasons. First, as the first analysis of planet sensitiv¬ 
ity from combined ground- and space-based microlensing ob¬ 
servations, it is natural to use the first single-lens event that 
allows us to do so. Second, Spitzer captured the peak of 


this event with a cadence higher than that of OGLE, and has 
an impact parameter (uq = 0.6) smaller than that of OGLE 
(no = 0.9), so even though Spitzer has fewer observations over¬ 
all, they are expected to contribute a significant fraction to 
the overall planet sensitivity. This helps us partly in testing 
the numerical result with our intuition based on microlensing 
theory, but mostly in understanding some related issues, as 
we will see below. Finally, the planet sensitivity of the first 
single-lens event with space-based parallax measurement has 
historical interest. 

We first construct four different planet sensitivity maps 
S(q,s) of OGLE-2014-BLG-0939 by treating this event un¬ 
der various different assumptions. First, we find S{q,s ) with 
parameters listed in Table 1, using both OGLE and Spitzer 
data, with no constraint; this sensitivity map is marked as 
“OGLE+Spitzer” in Figure 2. As Calchi Novati et al. (2015) 
has pointed out, one can find the constraint on the source 
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100 <x 2 <200 
200 <x 2 <500 
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Fig. 3.— The x 2 distributions of fake OGLE-2014-BLG-0939 light curves with planets placed at different positions (x,y). The black/red lines indicate the 
source trajectories seen by OGLE /Spitzer, with the dots representing the positions where the observations were taken. The lens is placed at (0,0). Planets in 
the upper panels have mass ratio q = 10 -2 , and planets in the lower panels have mass ratio q = 10 -3 . For both mass ratios we show the % 2 distributions from 
“OGLE+Spitzer,(—,—)” (the same as “OGLE+Spitzer,color” in Figure 2), “OGLE only”, and “OGLE+Spitzer,(+,—) cases. The grey solid circle in each panel 
is the Einstein ring, and the grey dashed circles have 1.35 and 1.55 Einstein radii, which correspond to the most sensitive s values for Spitzer data and OGLE 
data, respectively. The discontinuities in the \ 2 distribution in lower panels are caused by the discontinuity among observations. Note that the (+,—) solution 
(rightmost panels) is included here for pedagogical purpose only and is excluded by the measurement of the lens-source relative proper motion. We denote the 
injected planet as “detected” if the fake planetary light curve has x 2 > 200. 


OGLE-2014-B LG-0939 


log(a i /AU) 
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Fig. 4.— Comparsion of planet senstivities for two degenerate geometries that might in principle be inferred from the same OGLE-2014-BLG-0939 light curve. 
Coding of planet sensitivity curves is similar to Figure 2. The solid, dashed, and dot-dashed curves represent the sensitivity curves from “OGLE+Spitzer,(—)” 
(the same as “OGLE+Spitzer,color” in Figure 2), “OGLE+Spitzer,(+,—)”, and “OGLE only” (not shown in the left panel), respectively. In fact, the (+,-) solution 
was ruled out by Yee et al. (2015a), but is shown here to illustrate the impact of this type of degeneracy. 
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flux in Spitzer 3.6 //in band by doing a linear regression be¬ 
tween V — I and /- [3.6 /im] colors for non-microlensed field 
stars, and this flux constraint helped in a few cases to break 
the four-fold degeneracy. Therefore, we also compute S(q,s) 
with this flux constraint imposed, with the result marked as 
“OGLE+Spitzer,color” in Figure 2, and consider this Siq.s) 
as the final planet sensitivity map of this event and the stan¬ 
dard against which the others should be compared. In the third 
test, we use only OGLE data, and the S(q,s ) based on this is 
called “OGLE only”. 

In the above cases we treat this event as selected for 
Spitzer observations objectively. However, as has been 
shown in Yee et al. (2015b), events can also be selected for 
space-based parallax measurements subjectively. This is the 
case for OGLE-2014-BLG-0939, because there was no pre¬ 
determined objective criteria for 2014 season events. There¬ 
fore, we also consider the case in which this event was se¬ 
lected subjectively. We set the date when it was chosen for 
Spitzer observations to be the public announcement date, or 
the alert date (Yee et al. 2015b). For OGLE-2015-BLG-0939, 
faiert = 6811. We then compute S{q,s ) following the proce¬ 
dure proposed by Yee et al. (2015b), which we described in 
Section 2. For each a at the chosen (q,s), we first fit the 
fake data that were released before f a ] ert and find x| L alert . If 

Xsl alert > 10, we regard the injected planet as having been no¬ 
ticeable, and thus reject this a; otherwise we fit the whole data 
set to find yi?, and compare it with the threshold we choose 
to determine the detectability at this a. The sensitivity map 
computed in this way is marked as “OGLE+Spitzer,alert”. 

The most prominent feature in Figure 2 is the two separated 
triangular structures that are nearly symmetric about s = 1, 
which is a feature of low-magnification events (see, e.g.. Fig¬ 
ure 8 of Gaudi et al. 2002). The second remarkable feature 
in all but the “OGLE only” S(q,s ) maps is the double-peak 
structure at the bottom of each triangle diagram. This is better 
illustrated in Figure 3, in which we present the spatial dis¬ 
tributions of x 2 at various planet positions for two chosen q 
values. The reason for this feature is that the event was seen 
from two observatories to have fairly different impact param¬ 
eters (0.6 vs. 0.9) leading to two values of s at which the 
planet sensitivity peaks. For low and intermediate magnifica¬ 
tion events, the value of s at which the S(q,s) map peaks is set 
by (Gaudi et al. 2002) 


s — 


1 

s 


= M() . 


This gives us log.? = ±0.13 for Spitzer, and log.v = ±0.19 for 
OGLE, both consistent with our numerical result. 

As seen in panel (b) of Figure 2, the planet sensitivity 
changes dramatically from excluding to including the Spitzer 
light curve, as has been expected. We find no noticeable im¬ 
provement in planet sensitivity by imposing the color con¬ 
straint, and no significant loss of sensitivity when treating 
this event as subjectively selected. The reason for both is 
that the Spitzer light curve captured the peak of the event: 
Calchi Novati et al. (2015) showed that only in events in 
which Spitzer did not capture the peak could the color in¬ 
formation improve the single-lens fit. The same applies to 
the planet sensitivity. Because the peak contributes most of 
the planet sensitivity, excluding data before f a i ert < to^piner 
can only marginally modify the sensitivity curves by re¬ 
ducing sensitivity to planets farther from the lensing zone 
(0.6-1.6R E ). 


The four-fold degeneracy is a generic feature of single-lens 
events if observed by two well-separated observers (Refsdal 
1966; Gould 1994). Although it was broken in the present 
case by the measurement of the source proper motion, we con¬ 
sider here what the planet sensitivity would be if this four-fold 
degeneracy had persisted. Since this four-fold degeneracy re¬ 
duces to a two-fold degeneracy when one ignores the direction 
of motion and only consider the mass and distance of the lens 
system, one only need consider two of the four solutions. In 
the case of OGLE-2014-BLG-0939, these are the (—,—) solu¬ 
tion, i.e., the correct solution, and its counterpart (+,—) solu¬ 
tion. 9 The (-,-) solution has Earth and Spitzer trajectories 
pass from the same side of the lens, whereas the (+,-) solu¬ 
tion has the two trajectories pass from different sides of the 
lens. See Figure 3 for an illustration of the source trajecto¬ 
ries seen from Earth and Spitzer for these two different so¬ 
lutions. Figure 4 shows the sensitivity map from this (+,—) 
solution in comparison with that from the (—,—) solution (i.e., 
“OGLE+Spitzer,color” in Figure 2) and that of “OGLE only”. 
It shows that the planet sensitivity from the (+,-) solution is 
considerably higher for giant planets (q > 10 -3 ) with respect 
to the (-,-) solution, and almost doubles that from “OGLE 
only”. The reason is that the OGLE and Spitzer light curves 
are sensitive to different a values in the (+,-) solution. In 
cases for which the degeneracy cannot be resolved, we rec¬ 
ommend taking the mean of the sensitivities from these two 
solutions, weighted by the likelihood for each solution to be 
the correct one, for example the x 2 difference and the Rich 
argument (Calchi Novati et al. 2015). 

3.2. Planetary event OGLE-2014-BLG-0124 

OGLE-2014-BLG-0124 is the first microlensing plan¬ 
etary system with a space-based parallax measurement 
(Udalski et al. 2015). About 4 kpc away from Earth, the 
planet has mass 0.5 Mj, and is separated from a 0,7M n star 
by flJ _«3.1AU. 

The event OGLE-2014-BLG-0124 lies in OGLE field 
BLG512, meaning that it was observed at OGLE’s highest 
cadence, about once every 20 minutes. Again, we only in¬ 
clude data that were taken after HID' = 6000 but before 
the end of the 2014 season, and find in total 6647 OGLE 
observations. In addition, this event also falls into the mi¬ 
crolensing fields of the Microlensing Observations in Astro¬ 
physics (MOA, Sako et al. 2008) and the Wise Observatory 
(Shvartzvald & Maoz 2012) surveys, and therefore received 
8865 and 1010 observations from MOA and Wise respec¬ 
tively during that time interval. Since this is a planetary event, 
intensive follow-up observations were taken after the planet 
anomaly was noticed. 19 However, as has been emphasized 
before, these follow-up observations cannot be used for the 
planet sensitivity analysis. We mark all the ground-based ob¬ 
servations from the three survey groups as the “ground” data 
set. 

OGLE-2014-BLG-0124 received in total 39 Spitzer obser- 

9 A better way to mark the four solutions was introduced in Zhu et al. 
(2015a). However, this old label system is still used here in order to keep the 
consistency with Yeeetal. (2015a). The four solutions (—,+), (+,+)> (—. —), 
and (+,-) by the definition of Yeeetal. (2015a) and Calchi Novati et al. 
(2015) correspond to (+,+), f+, —), (—,—), and (—,+) solutions by the definition 
of Zhu et al. (2015a), respectively. 

10 Among all follow-up teams, only the Microlensing Follow-Up Network 
(/rFUN, Gould et al. 2010) obtained a few observations before this event be¬ 
came anomalous, but the purpose of these a few observations was to better 
characterize the source rather than to find planets. 
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Fig. 5.— Planet sensitivity curves of OGLE-2014-BLG-0124. The left panel shows the S(q,s ) map, and the right panel shows the S(q) distribution. The solid, 
dashed, and dotted curves represent the sensitivity curves from “OGLE+Spitzer”, “OGLE only”, and “OGLE+Spitzer,alert”, respectively. The position of the 
detected planet OGLE-2014-BLG-0124Lb is indicated as a blue dot with error bar in the left panel and a vertical blue dashed line in the right panel. The projected 
separation s = 1 in this case corresponds to a physical separation of 3.1 AU. 


vations during the interval HID' = 6815.5 to 6850.6. See Fig¬ 
ure 1 of Udalski et al. (2015) for the distribution of these ob¬ 
servations. Although these Spitzer observations play a signif¬ 
icant role in measuring the parallax effect, as Udalski et al. 
(2015) has shown, they are expected to contribute only a tiny 
fraction to the overall planet sensitivity because the trajec¬ 
tory of the source relative to the lens is very similar as seen 
from the ground and from Spitzer. Even so, however, we can¬ 
not include all 39 Spitzer observations for the planet sensi¬ 
tivity analysis, because some of them were carried out in re¬ 
sponse to the planetary anomaly found in the ground-based 
data. Considering that the Spitzer cadence was about once per 
day before the planet anomaly alert, we remove some Spitzer 
observations to form a time sequence with similar cadence. In 
the end, 29 Spitzer observations are included. Some example 
fake light curves generated from these two data sets are shown 
in Figure 1 . 

We construct three sensitivity maps for OGLE-2014-BLG- 
0124: one with only ground-based data, one with both ground 
and Spitzer data but no color constraint, and one by treating it 
as subjectively selected. We do not consider the impact of the 
color constraint, since the previous case has shown that if the 
peak of the event seen from Spitzer was captured, the color 
constraint can only modify the sensitivity curve marginally. 
We do not consider the case of four-fold degenerate solutions 
here either, because the opposite solution, which would have 
tte,n ~ 0.35, could be effectively ruled out by the ground- 
based data alone even if the planet were not detected. As 
with OGLE-2014-BLG-0939, in the case in which we treat 
the event as subjectively chosen, the date when this event was 
chosen for Spitzer observation is used as the alert date, i.e., 

falert = 68 11. 

Figure 5 shows these sensitivity maps S(q,s) and marginal¬ 
ized sensitivities S(q), with the position of the detected planet 
OGLE-2014-BLG-0124Lb (q = 7 x 10“ 3 and s = 0.94) marked 
for reference. As expected, although Spitzer observations are 
key to measuring the parallax effect, the removal of the entire 
Spitzer light curve has negligible effect on the overall planet 


sensitivity. In the case that this event was selected subjec¬ 
tively, the sensitivity drops by up to ~ 30%, but still remains 
significant, because the alert date 4i ert was fairly far from the 
event peaks as seen from Earth. 

4. DISCUSSION 

We present the first planet sensitivity analysis for mi- 
crolensing events with simultaneous observations taken from 
space and ground. We apply our methodology to two 
such events, the single-lens event OGLE-2014-BLG-0939 
(Yee et al. 2015a) and the planetary event OGLE-2014-BLG- 
0124 (Udalski et al. 2015). Assuming an underlying planet 
population of one planet per dex 2 per star, we find that OGLE- 
2014-BLG-0939 with a lens distance of 3.1 kpc shows an 
overall 11.5% detection efficiency for planets with mass ra¬ 
tio in the range 10 -4 < q < 10 -2 (or mass range from 6AW 
to 2Mj), and that OGLE-2014-BLG-0124 with a lens dis¬ 
tance of 4 kpc has 45% detection efficiency for planets with 
10 -5 < q < 10 -2 (or mass range from 2M® to 7Mj). 

The contributions to the overall planet sensitivity from 
space-based observations are considerably different in the two 
events. For OGLE-2014-BLG-0939, although Spitzer took 
fewer observations than OGLE did, the space and ground ob¬ 
servations contribute nearly equally to the overall sensitiv¬ 
ity, because the event had higher magnification as seen from 
Spitzer and was observed by Spitzer with a higher cadence. 
This result means that the planet sensitivity could have been 
higher if the solution with the source passing by the lens from 
different sides as seen from Spitzer and Earth were allowed, 
because then Spitzer would probe a significantly different re¬ 
gion surrounding the lens (see the rightmost panels of Fig¬ 
ure 3). It also suggests that a dedicated microlensing survey 
from space, when combined with ground-based surveys, can 
dramatically increase the probability to find planets, besides 
its capability to efficiently measure the microlens parallax ef¬ 
fect (Gould & Horne 2013). For OGLE-2014-BLG-0124, the 
ground-based observations are much more intensive and cover 
a much longer baseline, leading to a negligible contribution to 
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the planet sensitivity from Spitzer data, although these Spitzer 
observations are key to measuring the microlens parallax. 

Although neither of these two events showed high- 
magnification behavior, they both showed substantial planet 
sensitivities as seen in Figures 2 and 5. This result has 
implications for conducting microlensing experiments using 
narrow-angle space-based satellites such as Spitzer. Lim¬ 
ited by their field-of-view, such narrow-angle satellites can 
only be used to follow up events that were found from the 
ground. Ideally, one would want to follow up events that 
are very sensitive to planets, for example, high-magnification 
(A max > 100) events. However, two intrinsic difficulties of 
high-magnification events render such a strategy difficult to 
implement. First, high-magnification events are intrinsically 
very rare: one will only expect to have a few within the ~ 38 
day Bulge window for the Spitzer telescope. Second, due to 
operational constraints, using Spitzer -like telescopes for mi¬ 
crolens parallax measurements requires one to choose events 
at least a few days before the observations start. 11 However, 
it is difficult to confidently identify events that will achieve 
high magnification in advance. Our result here suggests that 
an ensemble of the more common low to moderate magnifi¬ 
cation events can also yield significant planet sensitivity and 
therefore probability to detect planets. 

Another key component of doing microlensing experiments 
using telescopes such as Spitzer is event selection. To make a 
well controlled sample of events suitable for statistical stud¬ 
ies, it is best to select all events based on some objective crite¬ 
ria. However, subjectively choosing events is also necessary 
for two good reasons. First, the objective criteria cannot cap¬ 
ture all events of interest, and second, an earlier subjective 
trigger may make the difference between measuring a paral¬ 
lax or not (Yee et al. 2015b). As illustrated in Figure 2c, the 
planet sensitivity of OGLE-2014-BLG-0939 considering only 
the data after the event is subjectively selected is comparable 
to the full planet sensitivity. In this case, the reason is that 
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a significant fraction of the planet sensitivity comes from the 
Spitzer data themselves. However, in general, this illustrates 
the value of subjective event selections even though only part 
of the light curve may be considered when calculating planet 
sensitivity. 

The analysis of these two events illustrates the role of 
space-based observations of microlensing events in measur¬ 
ing planet sensitivity, which is one component of measuring 
the Galactic distribution of planets. First, the additional ob¬ 
servations from well-separated observatory probe a different 
region of the lens system, increasing the planet sensitivity. 
Second, this demonstrates that low and moderate magnifica¬ 
tion events show substantial planet sensitivity. Both of these 
points apply to space-based microlensing parallax observa¬ 
tions in general including the K2 microlensing campaign. Fi¬ 
nally, we have shown that even though events may be selected 
subjectively for a targeted campaign like Spitzer s, they may 
still contribute significant planet sensitivity to the final sample 
as suggested by Yee et al. (2015b). 
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